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We investigate the transport of dipolar indirect excitons along the growth plane of polar 
(Al,Ga)N/GaN quantum well structures by means of spatially- and time-resolved photoluminescence 
spectroscopy. The transport in these strongly disordered quantum wells is activated by dipole-dipole 
repulsion. The latter induces an emission blue shift that increases linearly with exciton density, 
whereas the radiative recombination rate increases exponentially. Under continuous, localized ex¬ 
citation, we measure a continuous red shift of the emission, as excitons propagate away from the 
excitation spot. This shift corresponds to a steady-state gradient of exciton density, measured 
over several tens of micrometers. Time-resolved micro-photoluminescence experiments provide in¬ 
formation on the dynamics of recombination and transport of dipolar excitons. We account for 
the ensemble of experimental results by solving the nonlinear drift-diffusion equation. Quantitative 
analysis suggests that in such structures, exciton propagation on the scale of 10 to 20 microns is 
mainly driven by diffusion, rather than by drift, due to the strong disorder and the presence of 
nonradiative defects. Secondary exciton creation, most probably by the intense higher-energy lu¬ 
minescence, guided along the sample plane, is shown to contribute to the exciton emission pattern 
on the scale up to 100 microns. The exciton propagation length is strongly temperature dependent, 
the emission being quenched beyond a critical distance governed by nonradiative recombination. 


I. INTRODUCTION 

In the last years, indirect excitons (IXs) have attracted 
considerable attention because they constitute a model- 
system for studies of cold gases of dipolar bosons in a 
solid-state environment 1 . An IX is composed of an elec¬ 
tron and a hole, bound by Coulomb attraction but engi¬ 
neered in such a way that that a potential barrier sep¬ 
arates them along one direction in space. One of the 
most interesting and most studied realizations of IXs 
has been obtained in (Al,Ga)As/GaAs coupled quantum 
wells (CQWs) under externally applied electric field 2 . In 
such structures the electric field applied along the struc¬ 
ture axis tilts the conduction and valence band edges, so 
that the optical recombination takes place between an 
electron and a hole confined in different quantum wells 
(QWs), which form an IX. For typical electric fields of 
several volts per micrometer, red shift (quantum confined 
Stark effect) reaches tens of meVPEl The separation be¬ 
tween electron and hole within an IX increases drastically 
its radiative lifetime. Radiative lifetimes as large as t ens 
of microseconds have been reported in the literature!^ 
It has been shown, that if the thermalization time of IXs 
is shorter than their lifetime, the IXs can cool down to 
temperatures comparable or smaller than the degener¬ 
acy temperature, thus allowing for condensation into a 
coherent state.®®! 

Another important feature of IXs, is their intrinsic 
dipole moment, aligned perpendicular to the QW plane. 
These dipoles constitute an ideal driving force for putting 
IXs in motion, since an IX cloud rapidly e xpand s along 
the QW plane under its own repulsive forces! 3i ®2lDipole- 
dipole interaction also leads to the screening of the elec¬ 
tric field and induces a density-dependent blue shift of the 


exciton energjffH, at least below the Mott density above 
which excitonic effects are washed outP^ Blue shifts of 
IX energy up to several meVs and propagation lengths 
reaching a hundred of microns have been o bser ved in 
GaAs/(Al,Ga)As CQWs at low temperatures.® 13 

In-plane motion of IXs can also be controlled by ex¬ 
ternal voltages. This can be used to collect and manip¬ 
ulate IXs in trapsp®EhD latticed- 1 ^ and narrow channels™ 
which are similar to t he t raps and lattices used in the 
studies of cold atoms!- 19 20 This degree of control opens 
new possibilities for studying fundamental physical prop¬ 
erties of IXs.^ The efficient control of IXs by external 
voltages also paves the way to novel optoelectronic de¬ 
vices based on excitonic properties, as recently demon¬ 
strated!^®! However, the development of such devices 
operating at room temperature seems difficult, at least 
with GaAs/(Al,Ga)As heterostructures. This is mainly 
due to the relatively small exciton binding energy. In¬ 
deed,exciton binding energy is equal to 4 meV in bulk 
GaAs, which corresponds to a temperature of 45 K.® It 
is enhanced by quantum confinement in QWs but it is 
still reduced down to a few meV in CQWs, because of 
the on-axis separation of the electron-hole pair. 

Wide band-gap semiconductors are, from this stand¬ 
point, much better candidates: binding energies as high 
as 25 meV, for bulk GaN, or 60 meV for bulk ZnO— ; 
permit the preservation of excitonic properties up to 
room temperature. These compounds, in their natural 
wurtzite crystal form, exhibit another interesting fea- 
ture: spontaneous and piezoelectric polarization ^ !— 

So far, this was considered detrimental for optoelectronic 
applications but it can constitute a decisive advantage 
for creating of cold gases of dipolar excitons. Indeed, 
GaN/(Al,Ga)N or ZnO/(Zn,Mg)0 single QWs grown 
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FIG. 1. Calculated electronic band structures showing energy 
and wave functions of electrons and holes in a polar single 
quantum well (sample C) at (a) low and (b) high exciton 
density respectively and at T= 4 K. Eaigo.n is the energy gap 
in the barrier, Eo (Ex) is the exciton recombination energy in 
the zero (high) density limit, Ebs is the population-induced 
blue shift Ebs = Ex — Eo. 

along the (0001) crystal axis, naturally exhibit built-in 
electric fields of hundreds of kV/cm. This results from 
the difference of t otal p olarizations between the well and 
barrier materials !--^ For instance, in GaN/(Al,Ga)N 
QWs wider than 3 nm, the quantum-confined Stark effect 
is so strong that it dominates over the quantum confine¬ 
ment. In this case, the ground-state exciton energy lies 
below the band gap of GaN (3.48 eV at 4 K)P2 

In such QWs, IXs are naturally created in the ab¬ 
sence of an external electric field. Their binding energy 
decreases significantly when the well thickness increases 
but still remains largely above 10 meVP 1 ^ Consistently, 
radiative lifetimes increase exponentially and values as 
large as tens of microseconds h ave been measured for 
QWs widths of 6 — 8 nm only.EC^^fi] Q ne important con¬ 
sequence of the large band offsets, lifetimes and electric 
field effects in polar GaN/(Al,Ga)N QWs is the possi¬ 
bility to accumulate dipolar exciton densities larger than 
10 12 cnr 2 using reasonably small laser power densities, 
inducing blue shifts of several hundreds of meV. This sit¬ 
uation has been both observed experimentally and thor¬ 
oughly modeled,^ as shown in Fig. [lj By solving self- 
consistently the Poisson and Schrodinger equations, the 
curvatures of both potential and envelope wave-functions 
can be calculated, to obtain the optical transition ener¬ 
gies. 

The purpose of this paper is to exploit dipole-dipole re¬ 
pulsion to study the transport of IXs in GaN/(Al,Ga)N 
QWs grown along the (0001) axis. The tightly localized, 
strong nonresonant excitation in such structures pro¬ 
duces a potential that pushes the dipolar excitons away 
from the excitation spot. Time and spatially-resolved 


photoluminescence (PL) measurements allow us to ob¬ 
serve the emission/propagation of those excitons over 
tens of micrometers along the sample plane, at low tem¬ 
perature (T = 4 K). When the temperature is increased, 
the dipolar exciton propagation exhibits a drastic shrink¬ 
ing of extension, in the form of a horizon which corre¬ 
sponds to the distance at which the radiative recombi¬ 
nation rate becomes smaller than the nonradiative one. 
Fittings of the experimental results in the framework of 
the drift-diffusion model suggests the domination of dif¬ 
fusion over the drift induced by dipole-dipole repulsion. 

II. SAMPLES 


Sample name 

A 

B 

C 

QW width (nm) 

6.0 

6.0 

7.8 

A1 (%) 

20 

20 

15 

Inner barrier 

Alo. 2 Gao. 8 N 

AIN 

- 

F (kV/cm) 

900 

900 

670 

Trad (ps) 

22.5 

24.0 

1.6 

4 >0 (eV x cm -2 ) 

6.5 x 10 -14 

7.49 x 10 -14 

5.39 x 10 -14 

7 (cm -2 ) 

9.08 x 10 11 

11.6 x 10 11 

5.68 x 10 11 

Eo (eV) 

3.08 

3.13 

3.20 


TABLE I. Sample parameters: GaN QW width, A1 content 
in GaAIN barriers, inner barrier composition, built-in electric 
field F. The values of r° ad , 7 (Eo, 4>o) are used to calculate 
density-dependent radiative lifetimes T ra d (exciton energies 
Ex), according to Eq. [ 2 ] (Eq. [l}. 


The three samples under study were grown by molec¬ 
ular beam epitaxy (MBE) on c-plane oriented sapphire 
substrates followed by a 2 /Ltm thick GaN buffer layer. 
Their main characteristics are summarized in Table U 
For sample A, the active zone consists of a 30 nm thick 
A-lo. 20 Gao. 80 N barrier, a 4.4 nm (17 atomic monolayers) 
GaN QW followed by one monolayer of Alo. 20 Gao. 80 N, 
a 1.6 nm (6 atomic monolayers) GaN QW and a 30 nm 
thick Alo. 20 Gao. 80 N cap layer. This sample can alterna¬ 
tively be described as a 6 nm-wide GaN/Al 0 . 2 oGa 0 . 8 oN 
QW, in which an ultrathin barrier of Alo. 20 Gao. 80 N has 
been placed at the location where the envelope functions 
of the electron and the hole overlap most weakly. In such 
structures, the overlap is further decreased and the ra¬ 
diative lifetime is enhanced, but with a limited impact 
on the exciton binding energy. Sample B was fabricated 
following the same concept of a 6 nm-wide QW contain¬ 
ing a 1-monolayer barrier, but the latter was made of 

AIN, instead of Alo. 20 Gao. 80 N. For sample C, the active 
zone consists of a 30 nm thick Alo. 15 Gao. 85 N barrier, 
a 7.8 nm (30 atomic ML) GaN QW and 30 nm thick 

AIO. 15 Gao.ssN cap layer. Sample C was also studied in 
Ref. [ 35 ] Its low-temperature exciton lifetime was deter- 
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mined as r® ad = 1.65 ps. For samples A and B, using 
similar experimental conditions as in Ref. (35) and a laser 
repetition rate of 4 kHz, we measured characteristic de¬ 
cay times of 22.5 and 24 [is respectively. As for Sample 
C (Ref. (35|, the PL spectra of Samples A and B peak at 
much lower energies than the bulk GaN exciton. They 
also exhibit intense optical-phonon replica, whic h is a n- 
other consequence of the strong on-axis dipole.^ 37 -^ 


II.1. Excitonic energy and radiative lifetime versus 
excitonic density 



Exciton density, n (cm 2 ) 

FIG. 2. Calculated dependence of radiative lifetime (a) and 
transition energy (b) versus exciton density, for samples A, 
B and C found by solving Schrodinger-Poisson equations (tri¬ 
angles, circles and squares, respectively). Continuous lines 
show the approximations by exponential (lifetimes: a) and 
linear laws (energies: b). 

To describe the transport of indirect excitons quantita¬ 
tively, we need to establish a correspondence between the 
density of dipolar excitons and the blue shift of the op¬ 
tical transition. For this purpose, we solve numerically 
the coupled Schrodinger-Poisson equations by the self- 
consistent procedure described in Ref. [35] Above a Mott 
density of n ss 10 12 cm -2 , many body effects and band 
gap renormalization are known to influence the density 
induced blue shift, as explained in Ref. [35] introducing 
a simplified, phenomenological model for the additional 
terms to the field screening effects. Recently, a more 
thorough discussion and modeling of those many-body 
interactions has been published . 39 In the measurements 


reported below, we will be only weakly concerned by den¬ 
sities exceeding 10 12 cm -2 . Thus, we use the phenomeno¬ 
logical approach of Ref. [35] to establish the density de¬ 
pendence of the blue shift. A typical result is shown in 
Figure [lj for low (a) and high (b) exciton densities. The 
blue shift Ebs from the zero-density exciton transition 
E 0 can be observed at high density, as well as the in¬ 
creasing overlap of the electron and hole wave functions. 
Note also the strongly curved potential profile. The local 
electric field is not constant along the growth axis, mak¬ 
ing our calculation far more accurate than others which 
assume a constant electric field, rigidly altered by static 
charges distributed as Dirac functions at each interfaces. 
The values of the built-in electric field F at zero exciton 
density are given in Table [T] for each sample. 

Using this approach, we calculate the blue shift Ebs of 
the optical transition as a function of the exciton density 
n for each sample (Fig. [2jb)) . The dependence between 
the two quantities is nearly linear and the exciton energy 
can be approximated by: 

Ex = Eq + Ebs — E 0 + <j> 0 n , (1) 

where the self-induced potential <f>o is the fitting param¬ 
eter. Thus, from the experimentally determined PL blue 
shift we can deduce the IX density. 

The radiative lifetime T ra d is inversely proportional to 
the squared overlap integral of the electron and hole en¬ 
velope functions, and can be calculated as a function of 
the exciton density (Fig. [2](a)). The resulting density- 
dependent lifetime for all the samples is well described 
by an exponential 

Trad ^ T° ad exp(-n/'y), (2) 

where 7 is the fitting parameter. r 9 ad is measured by 
macro-PL measurements of the ultimate decay time con¬ 
stants, measured at low temperature and at long delays. 
The resulting parameters are reported in Table |T] for the 
three samples. 

III. EXPERIMENTAL SETUP 

For micro-PL experiments, the samples were placed 
in a cryostat and excited by a 266 nm continuous-wave 
laser (Crylas/FQCW 266 — 50) through a microscope ob¬ 
jective. This objective collects the PL signal in such a 
way that the image of the sample surface is reproduced 
at the entrance of the spectrometer, where it is filtered 
by a vertical entrance slit. The spectrometer is equipped 
with a grating with 150 gr/mm blazed at 390 nm. The 
detector is a cooled CCD camera with 1024 x 256 pixels. 
The pixel size is 26 pm, which corresponds to 760 nm on 
the sample surface. This value matches the smallest laser 
spot diameter that our objective can provide (800 nm), 
thus allowing for an overall spatial resolution of 0.8 pm. 
This setup provides a direct spatial imaging of the lu¬ 
minescence with a spectral resolution. Two-dimensional 
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FIG. 3. Color map of the PL intensity taken at T = 4 K on (a) sample A, (b) sample B, and (c) Sample C. The abscissa r 
corresponds to the distance from the excitation spot. The ordinate corresponds to the energy spectrum. The color map (d) 
shows the PL peak wavelength as a function of the position in the (x , y) plane of Sample B. The excitation spot is located at 
the center of the map. 


images, with spectral resolution, were obtained by scan¬ 
ning horizontally the lens which ensures the focusing onto 
the same entrance slit. The excitation power density was 
estimated as Pq = 2.3 MW/cm 2 . 


PL experiments combining time and spatial resolution 
are based on the same setup, except that the detector 
is a Hamamatsu streak-camera (model C10910 equipped 
with an S20 photocathode and UV coupling optics), en¬ 
hanced for ultraviolet detection. This camera can reach 
an ultimate time-resolution of 1 ps but, for the time- 
windows explored in this study, this resolution was rather 
of 15 ps. Pulsed laser excitation at 260 nm, through the 
same microscope objective, is provided by the frequency- 
tripled tunable radiation from a Titanium-Sapphire cav¬ 
ity, with typical pulse duration of 150 fs and tunable 
repetition rate. For the present ultra-slow decays, we 
set this rate to 4 kHz with an average power density 0.1 
MW / cm 2 . 


IV. EXPERIMENTAL RESULTS 
IV. 1. Low-temperature PL 

Figure [3](a-c) shows typical images of spatially and 
spectrally resolved PL obtained at T = 4 K for Samples 
A, B and C. The PL intensity is measured at different 
distances r from the excitation spot (horizontal axis) and 
at different photon energies (vertical axis), in the pres¬ 
ence of the strongly focused excitation at r = 0 at power 
density Pq. The dominant PL line, at higher energy, 
corresponds to the zero-phonon excitonic recombination. 
Optical phonon (LO) replicas of this PL line appear at 
lower energies, separated by multiples of 91 meV. As far 
as the spatial dependence is concerned, these replica fol¬ 
low the energy variation of the main zero-phonon line. 
In the following we will focus on the zero-phonon emis¬ 
sion line. Close to the excitation spot, the PL energy 
is blue shifted up to values close to the GaN excitonic 
gap EcaN = 3.48 eV. The maximal blue shift, observed 
at r ~ 0, is larger than 0.2 eV for Samples A and B, 
see Figure [3|(a-b) , and close to 0.15 eV for Sample C, 
see Figure pfc). However, the determination of the ex¬ 
act blue shift at r < 3 /iin is less accurate than for the 
rest of the sample surface, because the observed PL spec- 
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FIG. 4. Grayscale map of the PL intensity taken on sample A for different temperatures: (a) 4 K, (b) 20 K, (c) 50 K and 
(d) 80 K. The solution of the diffusion equation with constant diffusion coefficient D and recombination rate R such that 
\JD/R = 12 fi m is displayed in (a) as a red dotted line. 


tra are quite broad, owing to multi-excitonic effects, and 
other many-body interactions.^ At r < 20 /im, the in¬ 
dependence of the PL maximum has the shape of an ar¬ 
row, characterized by a decrease of the blue shift accom¬ 
panied with a decrease of intensity (note the logarithmic 
color-encoding of the intensity). This emission profile 
will be referred to as the arrow-shaped pattern. 

At longer distances (r > 20 pm), the PL energy is 
nearly constant, but the signal persists at surprisingly 
long distances. For Samples A and B, we observe exci- 
tonic emission up to 100 pm away from the generation 
spot. For Sample C, with the smallest initial energy shift, 
the excitonic emission is still observable up to several tens 
of pm away from the excitation spot. 

Figure [3](d) displays a color map of the QW peak 
emission wavelength as a function of the real-space co¬ 
ordinates (x,y) on the surface of Sample B. The above- 
mentioned central arrow-shaped region appears to be of 
circular symmetry, with a radius of 20 pm. Then the 
PL peak wavelength remains nearly constant over dis¬ 
tances as large as 100 pm. Note that this wavelength 
is slightly altered by the presence of linear microscopic 
defects (cracks), most certainly owing to local variation 
of barrier composition and/or strain field, near such de¬ 
fects. Careful examination of Figure [3ja-c) also empha¬ 
sizes small but reproducible fluctuations of both intensity 
and PL energy shift, when moving away from the exci¬ 


tation spot. We attribute these features to local fluctua¬ 
tions of barrier composition and/or strain, which modify 
the local accumulation of excitons and therefore their en¬ 
ergy. 


IV.2. Temperature dependence of continuous-wave 

PL 

Figure [4] shows the temperature dependence of the 
spatially-resolved PL spectra for Sample A, obtained at 
power density Po- One can observe the shrinkage of the 
spatial extension of the PL, when the temperature is in¬ 
creased. The PL signal falls below the detection thresh¬ 
old of the setup above a critical distance. This critical 
point can be referred to as an emission horizon, and it 
gets closer to the excitation spot as the temperature is 
increased. 


IV.3. Time-resolved experiments 

To get a deeper insight into the dynamics of exci- 
ton propagation, we performed low-temperature time- 
resolved PL experiments. By conjugating the temporal 
slit of the streak-camera to a given point on the surface 
sample, we are able to record the PL intensity collected 
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FIG. 5. Time evolution of the PL intensity Iq for sample A, 
taken at T = 4 K and at (a)4 pm, (b)6 pm, (c)16 pm and (d): 
60 pm away from the excitation spot. The grayscale map in¬ 
dicates the intensity of emission at a given energy (horizontal 
scale) and at a given time (vertical scale). Red curves show 
the calculated change of transition energy, with time, for each 
distance. 

at a given energy and at a given delay. Figure [5] shows ex¬ 
amples, in the case of Sample A, of such time-dependent 
PL spectra, gathered at 4 different distances from the 
excitation spot. The time-domain has been restricted, 
here, to the first 4 microseconds of decay, but we have 
also measured such decays over a full time range of 50 
ps, and we have examined, too, the first 500 ns after the 
excitation time. 

In the vicinity of the excitation spot (Figure [5j a)), the 
observed dynamics is the same as in spatially integrated 
experiments®® less than 50 ns after the excitation pulse, 
there is a sudden blue shift of the PL emission which re¬ 
sults from the fast creation of dipolar excitons. Then, 
the exciton population decays and the PL peak energy 
decreases, until it reaches a steady value. Only in the lat¬ 
ter regime, the intensity time-decay becomes exponential 
and provides us with the decay times that we gathered 
in Table H) As we increase the distance between the ob¬ 
servation and the excitation spots (Figure [5j)b) ), we ob¬ 
serve the same kind of behavior, but with a smaller initial 
blue shift: lower exciton densities are reached, consistent 
with the continuous-wave PL results. It can still be ar¬ 
gued that the time needed to propagate over 6 pm is 
shorter than the time resolution available in these ex¬ 
periments. However, the most puzzling results are those 
obtained above 15 pm (Figure [5](c-d)). Indeed, even at 
these large distances, the PL signal builds up within a 
few nanoseconds after laser excitation. This observation 
proves that, for distances larger than 15-20 pm, the ob¬ 
served exciton recombination can in no way arise from 
the propagation of those excitons, at least in the sense 
of a classical drift/diffusion process. As discussed in the 
next section, no reasonable diffusion coefficient can ex¬ 
plain such a fast arrival of excitons far away from the ex¬ 
citation spot. However, the quenching of the PL tail ob¬ 
served at high temperatures indicates the excitonic origin 


of this signal. We suggest that this observation results 
from the direct optical creation of secondary excitons, 
probably by the intense luminescence emitted under the 
laser spot and guided along the sample plane. 

V. DISCUSSION AND MODELING 
V.l. Origin of the arrow-shaped micro-PL pattern. 

The arrow-shaped pattern apparent in Fig. [3ja-c) is a 
signature of the decreasing density of excitons, when the 
distance from the excitation spot increases. The dipole- 
dipole repulsion between excitons pushes them away from 
the excitation spot. Exciton propagation, accompanied 
by recombination, results in a decreasing exciton den¬ 
sity. Generation, diffusion, drift, localization, radiative 
and nonradiative recombination are the competing mech¬ 
anisms that we need to unravel, in order to model quan¬ 
titatively our observations. We will assume that the blue 
shift directly reflects the local exciton density, whereas 
the PL intensity also results from the competition be¬ 
tween radiative and non-radiative recombination pro¬ 
cesses. The spatial dependence of the blue shift, observed 
less than 20 /im away from the excitation spot, appears 
to be temperature independent below T < 50 K, within 
our experimental resolution. Since drift mechanisms are 
typically temperature dependent because of the Einstein 
relation, we can first infer that either the thermalization 
is not complete below T =50 K for short propagation dis¬ 
tances, or that the exciton drift plays only a minor role 
here, due to the complex multiscale disorder potential. 

The simplest way to reproduce the arrow-shaped dif¬ 
fusion pattern is to consider the diffusion of excitons 
generated by a delta-peaked excitation pulse at r = 0. 
Let us assume a constant diffusion coefficient D and a 
constant recombination rate R. Far from the excitation 
spot, the steady state diffusion equation is reduced to 
DAn = Rn. This equation has an analytical solution: 
(j) 0 n = A x K a (ry/R/D), where A is a constant and K, 0 is 
the modified Bessel function o f the second kind. This so¬ 
lution with A = 0.1 eV and \JD/R ~ 12 /im is shown by 
the red curve in Fig.[4ja) that we superposed onto the PL 
map. Even though we know that this model is based on 
the un-satisfied assumption of a constant exciton lifetime 
and a constant diffusion coefficient, it reproduces reason¬ 
ably the overall shape of the PL spatial dependence. 

V.2. Origin of the PL tail 

A model based on purely diffusive propagation of IXs 
can, in principle, reproduce both the PL arrow-shape and 
the absence of temperature dependence, for r < 20 /im. 
On the other hand, the tail of the PL line, for r > 20 /iin 
cannot be reproduced by diffusion of excitons, because 
it would imply an anomalously high diffusion constant. 
Indeed, time-resolved measurements of PL at 60 /iin 
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from the excitation spot indicate that excitonic emission 
builds up in less than 100 ns after the excitation pulse 
(Fig. [5](d)). This cannot either be reproduced by the 
exciton drift, because of the weak exciton density gradi¬ 
ent observed at r > 20 /jm. We therefore suggest that 
a significant part of the higher-energy photons emitted 
around r = 0 by the QW itself is guided along the sam¬ 
ple plane and is capable of creating, at larger distances, 
secondary excitons. Let us now address the tempera¬ 
ture dependence of this emission and the reduction of the 
horizon. First of all, a common property of III-nitride- 
based heterostructures is a high density of nonradiative 
recombination centers (which are mainly threading dis¬ 
locations), of the order of 10 10 cm -2 . This is the case in 
the present heteroepitaxial QWs. Therefore, as the tem¬ 
perature delocalizes IXs, it also enhances the rate of their 
capture by those nonradiative traps. Given the very large 
values of the radiative lifetimes in our samples, at least 
in the unscreened regime at r > 20 fn n, the competition 
between slow radiative and faster nonradiative recombi¬ 
nations explains why PL signals are quenched, when the 
temperature is increased, for distances exceeding a cer¬ 
tain horizon. This horizon is simply determined by the 
exciton nonradiative recombination rate, which increases 
with the temperature. 


V.3. The transport equation 

In this section, we describe a formal model for the ex¬ 
citon propagation. In first approximation we neglect the 
thermalization effects and consider that excitons are close 
to the lattice temperature as soon as they leave the ex¬ 
citation spot. The equation for the in-plane transport of 
indirect excitons is®l 

f ) 77 

— = -V • J + G - Rn, (3) 

where n is the exciton density, G and R are the gener¬ 
ation and recombination rates, respectively, J is the IX 
current density, which can be split into drift and diffusion 
components: J = Jdrift + Jdiff- 

The diffusion current Jdiff is given by: 

Jdiff = —DVn (4) 

where D is the exciton diffusion coefficient. The random 
disorder potential, caused by alloy disorder in the barrier 
and by QW width fluctuations, hinders the in-plane ex¬ 
citon propagation along the GaN/AlGaN QW planed 
For high exciton densities, however, the potential fluc¬ 
tuations are screened and excitons, on average, become 
more mobile. This effect can be taken into account in 
the framework of a thermionic mo del which gives the 
following density dependence: 


-U 0 


(5) 


where Dq is the diffusion constant in the absence of dis¬ 
order and Uq is the amplitude of the disorder potential 
at zero density. This formulation implicitly accounts for 
the modification of the fluctuation potential by screening 
of the internal electric fieldWe estimate Uq ~ 30 meV, 
which corresponds to the energy shift induced by a fluc¬ 
tuation of one monolayer in the QW width. D 0 repre¬ 
sents the free exciton bulk diffusion coefficient, which is 
reached when the disorder potential Uq is negligible com¬ 
pared to the blue shift (f>on and/or the thermal energy 
k B T. 

The drift term Jdrift in Eq.[3]can be tentatively rewrit¬ 
ten as: 


Jdrift = -/z?zV(</> 0 n) (6) 


where the exciton mobility /z is connected to the diffusion 
coefficient D via the Einstein relation: 


M = 


D 

k^T' 


(7) 


In Eq. [6j the drift is governed by the self-induced poten¬ 
tial (j)Qn only. We expect this relation to be valid only 
when Uq <C 4>q n, i.e. when the disorder potential is effi¬ 
ciently screened. 

The generation term G which appears in Eq. [3] will be 
divided into two parts: G = G 0 + Gb g ■ The first term 
Go corresponds to the exciton generation under the laser 
spot: 


Gq = (2N p /RI) exp (~r 2 /R 2 0 ), 

where r is the distance from the center of the excita¬ 
tion spot, Rq= 500 nm is the spot radius, N p is the 
exciton generation rate. The second term Gb g is semi- 
phenomenological. We introduce it because any attempt 
to model the ultra-fast appearance of exciton recombi¬ 
nation at r > 20 /im within a drift-diffusion model us¬ 
ing reasonable diffusion coefficients is doomed to fail. 
Our current understanding of this phenomenon is that 
a significant part of the higher-energy photons emitted 
around r = 0 are most certainly guided along the sample 
plane, between the QW and the surface. Despite small 
absorption coefficients, some of those photons can be re¬ 
absorbed in the regions of r > 20 /zm and recycled, by 
thermalization, into low-energy excitons, with densities 
comparable to, or lower than 10 11 cm -2 . This explains 
the constancy of emission energy versus r, whereas the 
nearly constant intensity, for 20 < r < 100 /zm, is con¬ 
sistent with the weakness of the absorption coefficient, 
i.e. with the vanishing attenuation of the excitation. We 
therefore introduce a weak background source term Gb g 
to mimic this photon-recycling process. For simplicity, 
we take it independent of r. 

The recombination rate R can be expressed as a func¬ 
tion of the radiative and nonradiative recombination 
times T rad and T nrad : 

R = 1 1Trad T 1 /^~nrad- 


D = D 0 exp 


4>o n + k B T 


( 8 ) 
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FIG. 6. Grayscale map of the PL intensity calculated using Eq. [3] at different temperatures: 4, 20, 50 and 80 K as a function 
of the distance r from the excitation spot, and emission energy. The drift current given by Eq. [6] is taken into account. A 
phenomenological Gaussian broadening of the blue shift with a standard deviation of 70 meV has been introduced. 
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FIG. 7. Grayscale map of the PL intensity, calculated as in Fig. [6] except that only the diffusion current Jdiff (Eq. |4| is taken 
into account whereas the drift current Jdrift (Eq.§ is neglected. 


The density dependent radiative recombination time T ra d 
= T® a d ex P(~ n /l) is determined according to the pro¬ 
cedure described in the Section ED and does not de¬ 
pend on the temperature. By contrast, nonradiative time 
is temperature dependent. Indeed, when either exciton 
concentration or temperature increases, the excitons de- 
localization leads to more efficient capture of excitons on 
the nonradiative defects. These effect gives rise to the 
quenching of the PL at high temperature and determine 


the exciton nonradiative horizon. 

According to Refs. 03] and 0H the density of a exci¬ 
ton flux j flowing through a nonradiative recombination 
center of size a <C l, l the mean free path, is indepen¬ 
dent of the shape of the center and equal to j = nv/ 4, 
where v is the average absolute velocity of the excitons 
for which a uniform distribution of the velocity direction 
is assumed. Therefore, assuming that all excitons which 
reach the recombination center are captured, the rate of 
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exciton capture is proportional to 

1 l^nrad = — 27 TCl/l nrad ( 9 ) 

where l nra d is the average distance between the nonra- 
diative recombination centers. Using the relation D = 
v 2 t/2 , where r is the momentum relaxation time, we fi¬ 
nally obtain the relation: 

Tnrad = T° rad exp ^ bT ) ’ 

where T® rad , which we will use as an adjustable parame¬ 
ter, is the nonradiative recombination time in the absence 
of structural disorder Uq = 0. 

As the in-plane exciton propagation has a cylindrical 
symmetry centered at the excitation spot, Eq. [3] reduces 
to a one-dimensional partial differential equation. For the 
numerical resolution, we implemented a one-dimensional 
finite element scheme, where the non-linear Eq. [3| is re¬ 
placed by a linear one by using the Euler method. For 
the simulation of the time-resolved experiments, the tem¬ 
poral evolution of the exciton density is calculated within 
a time-implicit scheme, assuming generation rates Go and 
Gbg to be non-zero only at the initial time t = 0. For 
the experiments under constant illumination, we take 
dn/dt = 0 and G is time-independent. In both cases 
we use Neumann boundary conditions at r = 0 and 
r = 200 n m. Finally, the PL intensity is obtained from 
the exciton density as n/T ra d- 


Do 

T-0 

' nr ad 

U 0 

N p 

G bg 

1 

l 

30 

(a) 1.2 x 10 4 <5(i = 0) 

(b) 1.2 xlO 12 

N p x 10 8 


TABLE II. Main parameters of the fit for sample A. The 
diffusion constant is given in cm 2 /s, the recombination time in 
/rs, the disorder potential in meV. N p is the number of exciton 
generated per second for (a) time-resolved experiment and (b) 
experiments under continuous illumination. The background 
generation G bg is in m _ 2 s _1 and is taken proportional to N p . 
5 is the Dirac delta function. 


V.4. Numerical results and discussion 

The numerical implementation of the model points out 
that drift current is not the dominant transport mecha¬ 
nism for excitons. In all our attempts to introduce a drift 
current, increasing temperature resulted in a progressive 
shrinking of the spatial extension of the calculated arrow- 
shaped PL, and an increasing blue shift. Typical results 
of the simulation for Sample A are reported in Fig. [6] We 
have introduced the Gaussian broadening of the PL en¬ 
ergy with the standard deviation a = 70 meV, to account 
for the broadened PL lines. One can see that the arrow- 
shaped feature is significantly temperature-dependent, in 


contrast with the experimental result shown in Fig. [4] 
Similar results have been obtained for the two other sam¬ 
ples. 

The discrepancy with the experimental data suggests 
that Eq. [6] strongly overestimates the drift current, prob¬ 
ably because the gradient of exciton density induces too 
small energy variations compared to the disorder ampli¬ 
tude, which can present a complex multi-scale profile. 
By contrast, when only the diffusion term is included, 
the simulation gives results strikingly similar to the ex¬ 
periments. The solutions of Eq. [3j assuming Jdrift = 0, 
are presented in Figure [7] By using the set of fixed pa¬ 
rameters gathered in Table [I] and fitting parameters given 
in Table [XT] we reproduce the emission patterns observed 
experimentally. The onset of the nonradiative horizon, 
as the temperature is increased, is well accounted for by 
the thermionic description of the delocalization-capture 
of excitons by defects. 

This purely diffusive model also allows us to describe 
the time resolved experiments, using the same set of pa¬ 
rameters. The results of the calculation for Sample A 
is shown by red curves in Fig. [5] One can see that the 
onset of the PL blue shift associated which exciton diffu¬ 
sion from the excitation spot towards the detection point 
is probably too fast to be resolved in the present experi¬ 
ments (Fig.[5^a,b)). On the same time, the immediate (in 
less than 100 ns) onset of the excitonic emission at large 
distances from the excitation spot has the same origin 
as the PL tail in the spatially resolved images. It is due 
to secondary absorption of the PL guided in the sample 
plane, described by the semi-phenomenological genera¬ 
tion term G bg - 

The arrow-shaped patterns of PL measured and calcu¬ 
lated here are similar to the patterns that have b een re- 
ported in GaAs CQWs under external gate voltagd^SIll]. 
However, the typical energy and distance scale are quite 
different in our structures. In particular, the energy scale 
of the arrow-shaped feature reaches 200 meV, exceeding 
by more than an order of magnitude the GaAs values. 
Modeling of the PL patterns within the exciton drift- 
diffusion model taking into account the particularities 
of the nitrides, suggests that even in the highest den¬ 
sity and lowest temperature regime exciton transport is 
dominated by density-activated diffusion, rather than by 
the drift current. This effect seems to be proper to GaN 
QWs, but further studies are necessary for its better un¬ 
derstanding. In particular the design of heterostructures 
without guided PL and the resulting emission would be 
useful. 

It is instructive to analyze the density dependence of 
the diffusion coefficient (Eq. [5]) and recombination times 
(Eq-i Eq. [To] ), assumed in this model. These quantities 
are shown in Fig. [8] for Sample A. The density scale (up¬ 
per x-axis) is reported together with the blue shift scale 
(lower rr-axis). One can see that at low temperatures 
(T = 4 K) the diffusion coefficient is far from being con¬ 
stant in the range of the experimentally observed energy 
shifts. This results precisely from the dipole-dipole re- 
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fected by the exciton density, due to the more important 
contribution of the thermal energy to the exciton delo¬ 
calization. Even at low density the recombination time 
never reaches the radiative limit (which is temperature 
independent). This leads to the onset of the nonradiative 
horizon, well reproduced by this model. Finally, while 
both D and r change significantly with n, their product 
Dt changes much less. That is why the simple model of 
constant diffusion coefficient and constant recombination 
time presented in Fig. |4j(a) could successfully fit at least 
the PL blue shift. 


VI. CONCLUSIONS 


FIG. 8. Illustration of the modeling assumptions for Sam¬ 
ple A. Diffusion coefficient (a) and exciton lifetime r (b) as 
a function of the blue shift at T= 4 K (red solid line) and 
T= 80 K (red dashed line). Radiative and nonradiative com¬ 
ponents T r ad and Tnrad are shown by blue and green lines, 
respectively. 


pulsion energy, Ebs = 4>qti. The diffusion becomes expo¬ 
nentially slow at low density and only when the disorder 
is screened by the dipole-dipole repulsion the diffusion 
coefficient D reaches asymptotically the maximum value 
of D = D 0 = 1 cm 2 /s. This behavior is smeared out 
at higher temperatures, because the thermal energy con¬ 
tributes to the delocalization. For high exciton densities, 
comparable to, or larger than 10 12 cm -2 , the random po¬ 
tential fluctuations Uq are overcome even at T = 80 K, 
because Ebs Uo,kBT. This corresponds to the sit¬ 
uation near r = 0 in our PL experiments, and explains 
why the change of the temperature has very little influ¬ 
ence on the observations in this region. As the exciton 
density decreases from 10 12 to 10 11 cm -2 , the diffusion 
coefficient is reduced by two orders of magnitude, for 
T = 4 K, owing to the relocalization of carriers. This is 
what occurs, as the distance r is increased : the situation 
goes continuously from the domination of mobile, delo¬ 
calized excitons to frozen excitons, localized at random 
fluctuations. If the temperature is increased, as shown in 
Fig. |8j thermal excitation takes over (fc^T > Uo,Ebs) 
and therefore the diffusion persists even at lower exciton 
densities. 

Similarly, Fig. ^b) shows the calculated exciton life¬ 
time r and its radiative and nonradiative components. 
Despite a much smaller radiative lifetime at high densi¬ 
ties, the radiative efficiency, there, is not maximum: the 
delocalization of carriers favors their nonradiative cap¬ 
ture. On the other hand, at T = 4 K, reducing the 
density dramatically increases the radiative lifetime but 
increases the nonradiative lifetime even more, due to the 
relocalization of excitons. Therefore, at T = 4 K, the 
radiative efficiency is close to 1 for low exciton densities 
and the effective lifetime equals the radiative lifetime. 
At higher temperature the nonradiative time is less af- 


In conclusion, we have studied GaN/(Al,Ga)N c-plane 
quantum well structures using space and time-resolved 
PL spectroscopy. In the emission pattern that we repre¬ 
sent in the distance/energy coordinate space two distinct 
features are identified: (1) an arrow-shaped pattern up to 
20 (im from the excitation spot, which does not depend 
on the temperature (up to 80 K) and (2) a long range 
tail of the emission at constant energy, which disappears 
when the temperature is increased. The combination 
of spectral, spatial and temporal imaging of the exciton 
emission, and a detailed model of the transport and re¬ 
combination processes, allow to determine the most rele¬ 
vant phenomena explaining the short- and long-distance 
limits, and the short and long timescales. The arrow- 
shaped pattern is due to the density-activated exciton 
diffusion away from the excitation spot, and the density 
dependent emission energy. In contrast with GaAs-based 
structures, exciton transport is dominated by density- 
activated diffusion, rather than by the drift current. The 
long-range PL tail is attributed to secondary excitons, di¬ 
rectly created at long distances by higher-energy photons 
emitted around r = 0 by the QW itself and guided along 
the sample plane. This emission is efficiently quenched 
at high temperature beyond a certain distance, that we 
identify as the nonradiative horizon!- 24 

The results obtained in this work suggest that 
GaN/AlGaN QWs constitute a very promising system for 
the study of dipolar excitons and exciton-based devices. 
We have demonstrated that exciton transport takes place 
over more than 10 fim up to 80 K. Reducing the dislo¬ 
cation density by using GaN rather than sapphire sub¬ 
strate could substantially increase this distance. Conse¬ 
quently, such structures are promising for the studies of 
cold exciton gases and related quantum phenomena. The 
possibility to spatially separate the generation spot from 
the cold excitons reservoir could be crucial to achieve 
the quantum regime. In terms of the transport mech¬ 
anisms, it could be interesting to develop new models, 
taking into account energy relaxation and the formation 
of the phonon replica. The latter phenomenon is quite 
strong in the wurtzite semiconductors, compared to their 
cubic counterparts, and might affect the propagation of 
excitons. Other polar wide-band gap materials, such as 






ZnO/ZnMgO with excitons stable at room temperature 
is another possible research direction. 
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